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ABSTRACT 

Minor mergers of galaxies are expected to be common in a hierarchical cosmology such as ACDM. 
Though less disruptive than major mergers, minor mergers are more frequent and thus have the 
potential to affect galactic structure significantly. In this paper we dissect the case-by-case outcome 
from a set of numerical simulations of a single satellite elliptical galaxy accreting onto a massive 
elliptical galaxy. We take care to explore cosmologically relevant orbital parameters and to set up 
realistic initial galaxy models that include all three relevant dynamical components: dark matter 
halos, stellar bulges, and central massive black holes. The effects of several different parameters 
are considered, including orbital energy and angular momentum, satellite density and inner density 
profile, satellite-to-host mass ratio, and presence of a black hole at the center of the host. Black holes 
play a crucial role in protecting the shallow stellar cores of the hosts, as satellites merging onto a host 
with a central black hole are more strongly disrupted than those merging onto hosts without black 
holes. Orbital parameters play an important role in determining the degree of disruption: satellites 
on less bound or more eccentric orbits are more easily destroyed than those on more bound or more 
circular orbits as a result of an increased number of pericentric passages and greater cumulative 
effects of gravitational shocking and tidal stripping. In addition, satellites with densities typical of 
faint elliptical galaxies are disrupted relatively easily, while denser satellites can survive much better 
in the tidal field of the host. Over the range of parameters explored, we find that the accretion of a 
single satellite elliptical galaxy can result in a broad variety of changes, in both signs, in the surface 
brightness profile and color of the central part of an elliptical galaxy. Our results show that detailed 
properties of the stellar components of merging satellites can strongly affect the properties of the 
remnants. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

There now exists a strong body of evidence that lu- 
minous elliptical galaxies are, in many ways, a distinct 
group of objects when compared to lower-luminosity 
elliptical galaxies. This dichotomy is seen not only 
in thei r global structure - for example, isophotal 
shapes (|Kormendv fc Be"n der 1993) and degree of rota - 
tional support (|Davies et all 119831: iBender et all 11992ft 
- but also in their central properties: giant ellipti- 
cals (M v < —22) have shallow central surface bright- 
ness profiles (— dlnE/dlni? < 0.26), little rotation 
but velocity anisotropy, and boxy isophotal shapes, 
whereas spiral bulges and fainter ellipticals (M v > 
—20.5) exhibit central power-law surface brightness pro- 
files (— dlnE/dlni? > 0.5), signific ant rotational sup- 
port, and disky isophotal shapes (lLauer et al.1 119951: 
Faber et a l. 1997; Ravindra nath et al.l I2001h IL auer et al.1 
20051 : iFerrarese et al.1l2006l ). 

Different environments and formation histories are 
plausibly responsible for these two classes of galaxies. 
The fainter ellipticals are found in all environments in- 
cluding the field; their dense power-law centers may be a 
signatur e of dissipation in gas-rich mergers of d isk galax- 
ies (e.g.. iFaber et al.lll997t iGenzel et al.ll2001f ). The lu- 
minous ellipticals, on the other hand, are predominately 
found in galaxy clusters, some as brightest cluster galax- 
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ies and cD galaxies at or near the centers of clusters. 
Their shallow central luminosity profiles and low-density 
cores appear particularly difficult to form and preserve 
in theoretical models. 

Black holes provide an attractive model for the initial 
formation of the cores: mergers of two cuspy (p oc r~ 2 ) 
stellar bulges with central black holes result in the sink- 
ing of the black holes to the center of the merger rem- 
nant by dyn amical friction and the formation of a black 
hole binary (|Begelman et al.lll980f ). This process, and 
the subsequent orbital decay ("hardening") of the bi- 
nary, result in energy transfer from the black holes to 
the stellar system, and the combined effect is to reduce 
the i nitial cusp to a profile that is nearly cored in projec- 
tion fE bisuza kTet al.l 11991b iQuinlan fc Hernquist) 119971 : 
iMilosavlievic fc Merrittll200lD . 

How the shallow surface brightness profiles of lumi- 
nous ellipticals evolve and are preserved (for the most 
part) in subsequent mergers, however, is an open ques- 
tion and is the focus of this paper. The stellar mass 
within the core radius (defined to be the radius where 
the local logarithmic slope of the surface brightness is 
—0.5) is of order 1% or less of the total stellar mass; a 
small amount of mass deposit or removal can therefore 
easily alter the inner profiles. To maintain such shallow 
stellar profiles, massive elliptical galaxies, once formed, 
must somehow prevent additional stellar mass from being 
brought into their central regions through the frequent 
satellite accretion that is a hallmark of cold dark matter- 
based cosmological models. In particular, cosmological 
iV-body simulations that model only dark matter have 
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revealed a rich spectrum of dense dark matter subhalos, 
some of which can survive complete tidal disruption and 
orbit within their host halos. A logical expectation is 
that some of the massive core ellipticals have accreted at 
least one lower-mass galaxy with cuspier central density 
over their lifetimes. It is therefore intriguing that surface 
brightness cores in luminous ellipticals are ubiquitous. 

In this paper, we report the results from a set of merger 
simulations that explore the effects of a single satellite 
accreting onto a massive elliptical galaxy. Our simula- 
tions use realistic galaxy models composed of all three 
components that are dynamically relevant for elliptical 
galaxies: dark matter halos, stellar spheroids, and cen- 
tral supermassive black holes (section [5]) . As we show, 
each component plays an important role in the outcome 
of the accretion events and therefore should be included. 
The dark matter halos set the gravitational potential well 
on large scales and determine the orbital evolution in the 
early stages of the merger (before the stellar bulges in- 
teract strongly). The bulges dominate the gravitational 
potential on smaller (~ kpc) scales and thus affect the 
orbital evolution at late times in the merger. On sub-kpc 
scales, supermassive black holes can dominate the poten- 
tial well and thus contribute to the structural evolution 
of satellites, sometimes quite strongly. Even though the 
simulations are dissipationless, the three components ex- 
change energy gravitationally with one another during 
the merger, leading to observable effects on the stellar 
component that would not be fully captured if any com- 
ponent were neglected. 

An additional feature of our simulations is that we 
choose the orbital parameters based on distributions of 
orbits determined from large cosmological dark-matter- 
only simulations. By applying a semi-analytic pre- 
scription of dynamical friction and tidal stripping (sec- 
tion 2.3), we are able to start our simulations not at the 
virial radius of the host, where the cosmological distri- 
butions of halo merger orbits are determined, but at the 
smaller scale of four effective radii of the host's bulge. 

Our analysis focuses broadly on two aspects of satel- 
lite accretion: orbital decay and structural evolution of 
the satellite during the merger process (section [3} and 
observable signatures of infalling satellites on the host 
galaxies (section . For the former, we investigate the 
dependence of an accreting satellite's fate on the ini- 
tial orbital energy and angular momentum (section 3.1), 
the satellite's density (section 3.2), the presence of cen- 
tral black holes (section 3.3) and dark matter halo (sec- 
tion 3.4), and other factors such as the satellite-to- host 
mass ratio and the inner density profile of the satellite 
(section 3.5). For the observable consequences of satellite 
accretion, we investigate a number of different properties 
including surface brightness profiles (section 4.1), colors 
and color gradients (section 4.2), and mass deposit at 
large radii and kinematic structures (section 4.3). 

2. NUMERICAL SIMULATIONS 

2.1. Initial Galaxy Models 

The initial host galaxies in our main set of simula- 
tions (see Table 1) contain three components: dark mat- 
ter halo, stellar bulge, and a central supermassive black 
hole, all assumed to be in mutual equilibrium in the to- 
tal gravitational potential. We also perform comparison 
runs (a) without black holes to quantify the degree of 



TABLE 1 

Summary of Satellite Accretion Simulations 



Run 

(1) 



N P 

(2) 



Vt/Vr 

(3) 



0) 



7*.s 
(5) 



BH? 

(6) 



R e /kpc e/pc 
(7) (8) 



Three-component (dark matter+star+black hole) Runs 



HI 


1.1 x 10 6 





1.5 


1.0 


host 


1.39 


33 


Hlhi 


5.5 x 10 6 





1.5 


1.0 


host 


1.39 


20 


Hln 


1.1 x 10 6 





1.5 


1.0 


no 


1.39 


33 


His 


1 1 x 10^ 


n 
u 


x.o 


x.o 




r.oy 


■ ).) 


LI 


1.1 x 10 6 





0.75 


1.0 


host 


1.39 


33 


H2 


1.1 x 10 6 


0.5 


1.5 


1.0 


host 


1.39 


33 


H2n 


1.1 x 10 6 


0.5 


1.5 


1.0 


no 


1.39 


33 


H2d 


1.1 x 10 6 


0.5 


1.5 


1.0 


host 


1.07 


33 


H2dd 


1.1 x 10 6 


0.5 


1.5 


1.0 


host 


0.823 


33 


H2ddd 


1.1 x 10 6 


0.5 


1.5 


1.0 


host 


0.700 


33 


H3 


1.1 x 10 6 


1.0 


1.5 


1.0 


host 


1.39 


33 


H3s 


1.1 x 10 6 


1.0 


1.5 


1.5 


host 


1.39 


33 


M3 


1.1 x 10 6 


1.0 


1.1 


1.0 


host 


1.39 


33 


L3 


1.1 x 10 6 


1.0 


0.75 


1.0 


host 


1.39 


33 


L3d 


1.1 x 10 6 


1.0 


0.75 


1.0 


host 


1.07 


33 


H4 


1.1 x 10 6 


2.0 


1.5 


1.0 


host 


1.39 


33 


H4s 


1.1 x 10 6 


2.0 


1.5 


1.5 


host 


1.39 


33 


L4 


1.1 x 10 6 


2.0 


0.75 


1.0 


host 


1.39 


33 


Two-component (star+black hole) Comparison Runs 


Tl 


5 x 10 5 





1.0 


1.0 


host 


1.39 


20 


T2 a 


5 x 10 5 





1.0 


1.0 


host 


0.708 


20 


T3s a 


5 x 10 5 





1.0 


1.5 


host 


1.00 


20 


Tl 


5 x 10 5 





1.0 


1.0 


host 


0.700 


20 


T5 


5 x 10 5 





1.0 


1.0 


both 


0.700 


20 


T6 


5 x 10 5 





1.0 


1.0 


both 


0.700 


3.3 


T7 


5 x 10 5 


0.25 


1.0 


1.0 


host 


1.39 


20 


T8 


5 x 10 5 


0.5 


1.0 


1.0 


host 


1.39 


20 



Note. — Description of columns: 

(1) Name of run in order of increasing orbital angular momentum. 
H,M,L indicate high, medium, low orbital energy; n indicates no black 
holes; d indicates denser satellites; s indicates steeper inner density 
slope for stellar component of satellite 

(2) Total number of simulation particles in the host galaxy (before 
truncation). The total dark matter-to-stellar mass ratio is 10:1 for 
each galaxy (host and satellite) in the 3-component runs 

(3) Ratio of initial tangential to radial orbital velocity (vt/v r ) of 
satellite 

(4) Ratio of initial orbital velocity to local (host) circular velocity 

(5) Inner slope of the satellite's stellar density profile 

(6) Which of progenitors has a black hole 

(7) Effective radius R e of the satellite's stellar bulge (in kpc) 

(8) Spline force softening for stellar, dark matter, and black hole 
particles (in parsecs) 

a satellite-to-host mass ratio of 0.03 rather than 0.1 

satellite destruction due to the presence of the hole; (b) 
with two black holes, one at the center of each merg- 
ing galaxy to study the effects due to (unequal mass) 
black hole binaries; and (c) without a dark matter halo 
to quantify the differences in satellite orbit and fate be- 
tween the more realistic three-component runs and ear- 
lier work that modeled only stellar bulge and ignored 
dark matter. Luminous galaxies have non-negligible 
amou nts of dark matter on scale s of the half-mass radius 
(e.g- JPadmanabhan et a l. 2004), which can lead to non- 
negligible effects on merging satellites. The initial setup 
for each of the three components is discussed below. 

2.1.1. Dark matter halos 
The initial dark matter halos and stellar bulges in 



our simulations both have 
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iTremaine et al.lli994[ ): for component i (either dark mat- 
ter or stars), the density profile is given by 

«w = ( 3 -^(i)" 7 '( 1+ i)" <4 " 7 "' <» 

There are three parameters in this model: 7 controls the 
inner density slope, a is a characteristic radius, and M 
is the total mass of the given component. Additionally, 
we impose an outer density cut-off, mimicking the effects 
of tidal truncation. Details of this truncation and sev- 
eral other properties of these models are provided in the 
Appendix. 

We assign the dark m atter halos 7dm = 1 and match 
the inn er density of th i slHernquistj |l990) profile to that 
of the iNavarro et "all (|1997l NFW) prof ile having the 
same mass within its virial radius (see Springel et al. 
120051 for details of this procedure). Since the outer re- 
gions of our halos are truncated, the mass distribution 
on important scales is virtually identical to NFW pro- 
files. The mass and radius in our non-cosmological sim- 
ulations can be rescaled arbitrarily, but for better phys- 
ical insight, we choose to quote all the results in this 
paper for a fiducial (untruncated) dark matter halo mass 
and radius of (Mdm<, o)=(3.33 x 1O 12 M , 33 kpc) and 
(3.33 x lO n Af , 13.33 kpc) for the host and satellite, re- 
spectively We generally consider 1:10 mergers in this pa- 
per since this mass ratio is small enough be in the minor 
merger regime but large enough for mergers to occur on 
reasonably short ti mescale (see Sec 2.3 below and, e.g., 
iTaffoni et al.1 12003 for calculations of merger timescales 
for a variety of mass ratios). We also investigate 1:33 
mergers to check the robustness of our predictions (see 
Sec 3.5). 

2.1.2. Stellar bulges 

The stellar bulges are self-consistently embedded in 
dark matter halos with mass Mdm = 10 M* before trun- 
cation. For the host galaxies, the initial stellar bulges 
are all given 7*^ = 0.5 to model the nearly flat central 
surface brightness profiles of massive ellipticals (dashed 
curves in Fig. [1]). To assign the scale radius a in equa- 
tion ([TJ) for a given stellar mass M*, we use the mean 
size-mass relation , R e oc M° ,e , for early-type galaxies 
(|Shen et al . 2003) determined from the Sloan Digital Sky 
Survey (S PSS: lYork et al.|[2000]) an d the stellar mass es- 
timates of Kauffm ann et al.l (|2003f ). Note that the re- 
lation between the effective radius R e and scale radius 
a depends on 7: R e /a = (2.358, 1.815, 1.276,0.744) for 
7 = (0.5, 1.0, 1.5, 2.0). For our fiducial mass scale, the 
stellar mass of the host is = 3.3 x 10 11 M , which 
corresponds to a gala xy of magnitude M v w —22 (e.g., 
iGebhardt et al.ll2003l ). For 7*^ = 0.5, the corresponding 
effective radius is R e — 5.5 kpc (a = 2.33 kpc), and the 
core radius is approximately 400 pc, near the upper end 
of the obser ved distribution of core sizes at our chosen 
stellar mass ([Lauer et al.ll2006l Fig. 5). 

For the satellite galaxies, the initial stellar bulges are 
given 7*3 = 1.0, resulting in initial surface brightness 
profiles (solid curves in Fig. [T]) that closely approximate 
the r 1 / 4 law (jde Vaucouleursl[l948f ) 3 . We also perform 
comparison simulations with 7* iS = 1.5 satellites (keep- 
ing R e fixed) and find increasing 7* jS to be somewhat 

3 The Hernquist (7 = 1) profile actually deviates from 



degenerate with making the satellites more compact by 
reducing a (see Sec. I3.5|) . Since tidal effects depend on 
the local ratio of the satellite and host densities, we find 
it important to take into account the scatter in the SDSS 
R e — M* relation, whic h has a dispersion of <7i n R ~ 0.26 
for early- type galaxies (Shen et al. 2003). For satellites 
of a given M*, we therefore consider galaxies on the SDSS 
R e — M+ relation as well as those with R e that are la and 
2a smaller, referring to them as "standard" , "lcr dense" , 
and "2(7 dense" satellites in subsequent discussions. 

Fig- Q] shows the stellar density (left panel) and surface 
brightness (right panel) profiles of our initial galaxy mod- 
els for the host (dashed curves), the 7* !S = 1.0 satellites 
(solid curves), and the 7*^ = 1.5 satellite (dot-dot-dot- 
dashed curve). The satellites are denser than the host 
for r < 300 pc, while the surface brightness of the satel- 
lites exceed that of the host galaxy for R < 150 — 300 
pc (and only by a factor of 3 at 10 pc for the standard 
satellite) . Fig. Q] also shows the stellar density profile for 
a galaxy with the same stellar mass as the host, assum- 
ing 7*^ = 1 an d R e following the SDSS relation (dotted 
line) . Its inner density is nearly identical to that of the 
satellite, which underlines the distinction between core 
and power-law galaxies at fixed luminosity. 

The line-of-sight stellar velocity dispersions are 240 km 
s _1 for the host and 140 km s _1 for the ( standa rd) satel- 
lite, placing the m on the Faber-Jackson (Il976f) and fun- 
dame ntal plane dDressler et al.lll987tlDiorgovski fc Davis! 
Il987| ) relations. 

2.1.3. Central black holes 

We model the black hole as a point mass 
with Msr/M* = 2 x 10~ 3 , sim ilar to the re- 
lation observed in galactic n u clei flM agorria n et al.l 
19981 iFerrarese fc MerritH 120001 : IGebhardt et all 120001 : 
Haring fc Rixi I2004D . The black hole exerts a dy- 
namical influence over the stars (and dark matter) 
within its sphere of influence, tbh, which we de- 
fine as the radius enclosing a mass in stars equal to 
twice that of the black hole: M+(< 7'bh) = 2A/bh- 
With this definition, the standard tbh = GMbh/c 2 
holds true for the isothermal sphere. For refer- 
ence, r BH /a = (0.123, 0.0675, 0.0258, 0.00402) for 7 = 
(0.5, 1.0, 1.5, 2.0), assuming M BH /M* = 2 x 10~ 3 . Note 
this scale can be much larger than the « pc scale often 
attributed to black holes: for our choice for host galaxy 
of a = 2.33 kpc, 7*, h = 0.5, and M* = 3.3 x 1O U M , 
we find tbh = 290 pc, comparable to the size of the sur- 
face brightness core in our galaxy model (dashed line in 
the right panel of Fig. [l}. In comparison, a profile with 
7*.h = 2 but with the same M* and R e as our 7*^ = 0.5 
profile has tbh = 7.7 pc, showing how sensitive the black 
hole's sphere of influence is on the stellar density profile 
under consideration. 

For cases with two black holes, another relevant scale is 
'"hard = G/i/4<7 2 , which is the separation at which the two 
black holes form a hard binary. This is the radius where 
the gravitational binding energy of the binary exceeds the 
kinetic energy of the stellar background. Note that for an 
equal-mass binary in an isothermal stellar background, 

the r 1 / 4 law at both small and large radii (Hcrnauist 1990; 
IBovlan -Kolchin et al. 200^) but is a good match to the r 1 / 4 profile 
for the region containing most of the mass. 
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Fig. 1. — Initial galaxy models for the simulations. Left: stellar mass density profiles of 7*, s = 1 (set of solid curves) and 7* iS = 1.5 
(dot-dot-dot-dashed curve) satellites. Also shown is the initial host stellar density profile (dashed curve) as well as the host profile predicted 
by the SDSS R e — M* relation if the host had a Hernquist density profile (dotted curve). Right: surface mass density profiles corresponding 
to the curves on the left. 

that the characteristic scale of the wandering, r w w 
vbh I '(Gp) 1 / 2 , can be somewhat large for a low-density 
core. We do observe wandering in our simulations, but 
the phase space structure is not affected in any important 
way above our force softening scale. 

2.3. Choice of orbital parameters 

In order to make meaningful predictions about satel- 
lite accretion, we would like to use realistic orbital prop- 
erties drawn from merger events in cosmological simu- 
lations. This task is non-trivial, however, since merg- 
ers of dark matter halos and mergers of galactic bulges 
are not the same thing: analyses of cosmological sim- 
ulations generally compute merger orbital properties at 
a galaxy's dark matter vir ial radius (e.g., iBensonl [20051 : 
iKhochfar fc BurkertJ |2006| ) , whereas we are mainly in- 
terested in the detailed dynamics among dark matter, 
stars, and black holes after the two merging galaxies' 
stellar bulges begin to interact. It has long been known 
that galactic cannibalism may b e important for the 
grow t h of central cluster galaxies jOstriker fe Tremaind 

lLauerl [~ 



Hiard = »"bh/8, and this is always an upper limit for the 
cases we are considering. 

2.2. Methodology and Numerical Effects 

All simulations presented ar e performed us ing the par- 
allelized tree code GADGET-2 (|Springelll2005h in pure N- 
body mode. The initial positions for the dark matter and 
star particles are drawn from their respective spherically 
symmetric density profile described in Sec. 2.1. The ini- 
tial velocities for each component are drawn from that 
component's distribution function, assuming the system 
to be spherical and isotropic (see iBovlan-Kolchin et al.l 
[20051 for details). In using this procedure, we make no 
simplifying assumptions about the velocity structure of 
the model components and the resulting initial condi- 
tions are initially in an equilibrium that is modified only 
by numerical effects. An example of the numerical sta- 
bility of our galaxy models is provided in the Appendix. 

We have performed resolution tests to determine if the 
results reported in Sec. 3 and 4 are sensitive to the num- 
ber of particles, force resolution, or force softening used 
in the simulations. One such test (run Hlhi in Table [l} 
uses 5 times more particles for both the bulge and halo 
and a force softening that is a factor 5 1 / 3 smaller than 
that of the standard runs. No differences are found be- 
tween this run and the equivalent run (HI) at our stan- 
dard resolution. All runs have equal particle masses for 
the stellar and dark matter particles to avoid numerical 
relaxation; accordingly, the force softening is equal for 
both particle species in all runs. We use the same force 
softening for the black hole as well. In all production 
runs, energy is conserved to better than 0.2%. 

For simulations with central black holes, gravitational 
brownian motion can be significant in our isolated galaxy 
models. The expectation is that the black hole and 
star particles will scatter off of each other, redistributing 
the particle energies toward equipartition: (vbh/c *) 2 = 
3JV- 1 M*/M B h (though see iChatteriee et all [20021 for a 
more complete discussion) . This velocity scale for reason- 
able particle numbers is not insignificant: for Mbh/-^* = 
2 x 10" 3 and iV* = 10 5 , then v B h/<t* = 0.122. Note 



Il975t lHausman fc Ostrikerl fl978l : 



1 119881 ) but also 



must be relatively ineffi cient due to t he lo ng merging 
timescales for satellites ([Merrittl [1983. 119851) . As a re- 
sult, there is a selection effect favoring satellites with 
higher masses that have correspondingly shorter dynam- 
ical friction timescales. For the same reason, accretion 
events that lead to merging with a central galaxy should 
be skewed towards lower angular momentum orbits than 
the average orbit for all halos accreted at i? v ir- 

We choose to start the simulations when the center of 
the satellite is at a distance of 4i? e of the host, a dis- 
tance considerably smaller than the virial radius i? v ir of 
the host. This starting point has the advantage that we 
are still in the regime in which the stellar bulges have not 
yet interacted strongly, yet it avoids the often-times long 
inspiral from i? v i r to 4 R e that can be reasonably handled 
analytically. To map the orbital parameters from i? v j r to 
4 R e , we have performed a series of semi-analytic calcula- 
tions in which we integrate numerically a satellite's orbit 
in the background potential of a host, taking into account 
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the processes of dynamical friction (causing orbital de- 
cay) and tidal stripping (removing mass from outside of 
the Roche limit). The initial orbital parameters at i? v ir 
are chosen to be typ ical of those measured from cosmo- 
logical simulations (|Bensonl 120051 : iKhochfar fc Burkertl 
2006): we vary around a "most probable" orbit, which 
we take to have v r » w c ; rc and vt ~ 0.7u c i rc . For decel- 
eration due to dynamical friction, we use the standard 
formula 

-^-Worb = -47TG 2 In A M sat Phost (< «orb) , (2) 

where phost(< Vorb) is the density of host halo par- 
ticles with velocities less than the bulk velocity of 
the satellite and In A is t he usual Coulomb logarithm 
dBinnev fc Tremaind Il987t ). [Following the results of 
iTavlor fc Babull (|200lD . we use A = 1 + M host /M sat (t).] 
The drag force is thus proportional to M s 2 at . Although 
many of the assumptions that go into deriving the Chan- 
drasekhar formula are not appropriate for simulations 
such as ours, previous work has shown that this equation 
provides a reasonable approximation for the orbital de- 
cay of satellites when coupled with expressions for tidal 
stripping and shocking (e.g.. IVelazquez fc White! 119991 ; 
ITavlor fc Baf^l200~fl : [Taffoni et al.ll2003f >. 

For tidal stripping, the relevant scale is the satellite's 
tidal radius i?tid, which is determined by the balance 
between the satellite's self-gravity and host's tides. At 
separation r from the host center, i? t id is approximately 
given by p sa t(-Rtid) ~ Phost( r ) an d the mass outside of 
i?tid tends to be stripped away on a local dynamical time. 
We model this process by integrating the satellite orbit 
with timestcps that are a fraction / of the instantaneous 
orbital period, calculating the tidal radius at each point, 
and removing a fraction / of the mass outside the tidal 
radius at each timestep. 

Besides the processes of dynamical friction and tidal 
stripping discussed above, impulsive heating from fast 
encounters is expected to lead to an additional ef- 
fect of gravitational shocking (lOstriker et al.1 Il972t 
iGnedin fc Ostriker! Il99l iGnedin et all 11999ft . This ef- 
fect occurs when the tidal field of the host on the satel- 
lite varies rapidly, e.g., when a globular cluster passes 
through a galactic disk or, in our case, when a satellite 
galaxy approaches the stellar bulge of the host. After ini- 
tially compressing the satellite, the shock heats the satel- 
lite and causes it to expand, rendering it more susceptible 
to disruption on further pericentric passages or by tidal 
stripping. The strength of the shock (assumed to have 
characteristic timescale r s hock) depends on the satellite's 
internal structure and dynamics: stars on orbits with 
angular frequencies uj <C T s hock tend to receive a boost 
in energy, while orbits with u tT_], are not strongly 
affected (|Weinberg|ll994D . Since our semi-analytic calcu- 
lation is to map orbital parameters between i? v j r and the 
first crossing of 4_R e , tidal shocking should be unimpor- 
tant and is not included. Our calculati on thus amounts 
to a s im plified version of the models of ITavlor fc Babul 
(|2001ft or lTaffoni et al.l(|2003ft . The effect of gravitational 
shocking, however, is seen in our numerical simulations, 
as discussed in Sec 3. 

The results of our semi-analytic calculations suggest 
that when the satellite-to-host mass ratio is small (< 0.1) 
and the initial orbital velocity at i? v ir is high [w or b > 



^circC-Rvir), where v^ ilc (r) = GM(< r)/r], only satellites 
on relatively radial orbits (with tangential-to-radial ve- 
locity ratio of Vt/v r < 1) decay quickly enough to be 
relevant, i.e., reaching AR e within half a Hubble time. 
Much of the parameter space is therefore eliminated. For 
the orbits that do decay rapidly enough, we find that the 
satellite velocity at 4 R e relative to the local circular ve- 
locity v C i rc at 4i? e typically ranges from v OI h/v C i rc = 0.7 
to 1.75. The specific angular momentum at 4i? e tends 
to be high: even though we consider only orbits with 
v t /v r < 1 at the virial radius, it is not uncommon for the 
tangential velocity to be greater than the radial velocity 
at 4i? e . As a result, we choose initial orbital energies 
ranging between forb /^circ 

= 0.75 and 1.5, and angular 
momenta ranging between Vt/v r = and 2.0 at AR e . 
The full set of initial conditions for our simulations is 
listed in Table [T] 

3. RESULTS: FATE OF SATELLITES 

In this section we investigate how the survival of a 
merging satellite depends on its orbit and density, as well 
as on the presence of a black hole and dark matter halo. 

3.1. Dependence on orbital angular momentum and 

energy 

The three panels in Fig. [2] compare the time evolution 
of the satellite orbital decay and mass loss for three sim- 
ulations that differ only in the initial orbital angular mo- 
mentum, where the initial tangential-to-radial velocity 
ratio ranges from Vt/v r — (purely radial) to v t /v r = 2 
(runs HI, H2, and H4 in Table 1). The orbital energy 
is fixed at forb/^circ = 1-5. The vertical axis plots the 
center-of-mass separation between the satellite and host 
(solid line) and the corresponding evolution of the mass 
within the satellite's initial R e (dashed line) and initial 
scale radius a (the quarter-mass radius for a Hernquist 
density profile; dotted line). 

Satellite mass losses are seen to follow distinct steps, 
with most changes due to strong tidal shocking followed 
by tidal stripping as the satellite passes through the peri- 
center of its orbit. The satellite mass stays relatively con- 
stant between pericentric approaches for all orbits. This 
is not unexpected, since both tidal stripping and gravita- 
tional shocking are most effective at the pericenter of an 
orbit. Fig. [2] shows that the tidal effects at the first two 
pericentric passages are significantly stronger for more 
radial orbits, leading to faster mass loss in the satellites. 
It also shows that the orbital decay of the satellite with 
low angular momentum subjects it to increasingly strong 
tides and large mass loss. The end result is that the ra- 
tio of the final to initial satellite stellar mass within its 
initial effective radius R e (dashed line) is only w 10% for 
the purely radial (vt/v r = 0) orbit and increases to 18% 
and 27% for v t /v r — 0.5 and 2, respectively. 

A prominent feature in Fig. [2] is the spikes in the 
M(< r, t) curves at each pericentric passage, which are 
caused by tidal compressio n from the host pot ential ac- 
companying a tidal shock. De kel et al.l (|2003l ) have in- 
vestigated this effect for a variety of profiles and find 
that satellites in a host potential with 7 < 1 experi- 
ences compression in the radial direction (i.e. along the 
line separating the centers of mass of the host and satel- 
lite) in addition to the ubiquitous compression in the 
perpendicular direction. For 7 > 1, the tidal interac- 
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Fig. 2. — Evolution of the center-of-mass separation between the satellite and host (solid) for three simulations that differ only in the 
initial orbital angular momentum; the ratio of tangential-to-radial velocities, v t /v r , increases from left to right. The initial orbital energy 
is fixed at v or -^ /v c i IC = 1.5. Also plotted is the evolution of the satellite's stellar mass M sa t(< r, t) interior to the satellite's initial effective 
radius R e (dashed) and scale radius a (dotted), normalized to the initial stellar mass within that radius. (The scale corresponding to the 
dashed and dotted curves is given on the right vertical axis.) From these plots it is apparent that most of the heating of the satellite bulge 
occurs during pericentric passages. Furthermore, there is significant tidal compression during each pericentric passage. Satellites on more 
radial orbits experience larger mass losses and fall to the center of the remnant faster. 
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Fig. 3. — Same as Fig. [2]but for three simulations that differ only in the initial orbital energy. The initial orbital angular momentum 
is fixed at vt = v r . More energetic orbits lead to larger mass losses; it also takes significantly longer time for satellites on these orbits to 
decay to the center of the remnant. 



tion tends to stretch the satellite in the radial direction 
while compressing it in the perpendicular direction. The 
host galaxies in our simulations have stellar bulges with 
7*.h = 0.5, so there should be a net compression of the 
satellite as it passes near the center of the host, which is 
indeed the case in Fig. O 

Fig. [3] is similar to Fig. [5] but compares results from 
three simulations of different orbital energies at fixed 
vt/v r = 1. The corresponding runs in Table 1 are L3, 
M3, and H3. The figure shows that satellites with larger 
initial velocity i> or b (i.e. more energetic and less bound 
orbits) merge more slowly, experience more pericentric 
passages, and suffer more severe mass loss overall. The 
amount of mass loss at each pericentric approach, how- 



ever, is relatively similar for different orbital energies, 
indicating that the primary factors for determining the 
destructive effect of gravitational shocks and tidal strip- 
ping on a given satellite are the number of pericentric 
passages and the distance of each pericentric passage. 

3.2. Dependence on satellite density 

Fig. [4] compares the orbital decay and mass loss of the 
satellite for three simulations that differ in the initial 
satellite density but are otherwise identical. The satel- 
lites in the three panels are initially assigned the three 
densities shown in Fig.[T]and discussed in Sec. 2.1.2. The 
sensitive dependence of the amount of satellite disruption 
on its initial density is clearly seen: the denser the satel- 
lite, the greater its degree of survival. The dashed lines 
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Fig. 4. — Similar to Figs. [2] and [3] but for three simulations that differ only in the initial stellar density of the satellite. Left: local 
R e — relation from SDSS early-type galaxies. Center: 1<t denser. Right: 2a denser. The runs are otherwise identical (vt/v r = 0.5, 
"orb/'f'ch-c = li with dark matter and black hole). The profiles are computed at apocenter after each pericentric passage. More compact 
and denser satellites are noticeably better preserved in both the inner and outer structure. 



show the evolution of the mass within the satellite's ini- 
tial R e relative to the initial mass within the same radius. 
The final value of this ratio changes from ~ 20% for the 
run with the standard density to ~ 45% for the la dense 
run to ~ 60% for the 2a dense run, showing that satellite 
density significantly affects the satellite survival on the 
scale of R e . A similar effect is evident for the evolution of 
mass within the satellite's scale radius a (dotted lines). 

The density of the satellite also influences the region 
from which the satellite's mass is stripped. The stan- 
dard density satellite ends up losing fractionally more 
mass within a than within i? e , while the 1 and 2a dense 
satellite lose fractionally more mass within R e than a. 
In other words, the central region of the standard den- 
sity satellite is more affected than the half-mass radius, 
presumably because of the black hole's influence. The 
compact satellites are more resilient and their outer re- 
gions are therefore relatively easier to strip than the inner 
portion. 

For a closer look at the changes in the internal struc- 
tures of the satellites as they sink towards the center of 
the host galaxy, we show in Fig. [5] the satellite's stel- 
lar density profiles p sa ,t(r) (upper) and circular velocity 
profiles UcircM (lower) after each pericentric passage for 
the same runs as in Fig. [4j The profiles are computed 
at apocenter after the correspondingly labeled pericen- 
tric passage, e.g., the curves labeled "first" correspond to 
profiles at apocenter after the first pericentric pass. The 
circular velocity profiles clearly illustrate that the tidal 
effects during the first pericentric passage are confined 
to a satellite's outer regions, while the later passages are 
responsible for the inner mass loss. The satellite with 
the standard density is stripped progressively more with 
each pericentric passage, resulting in a final central den- 
sity that is reduced by about two orders of magnitude 
and a maximum circular velocity u ma x that is reduced 
by 40% relative to the initial condition. The la dense 
satellite is much less affected in density (although v max 
is reduced by 30%), while the most compact (2a denser) 
satellite is stripped the least. 



The effect of a black hole in the host on the satellite 
can also be seen clearly in the two upper-right panels of 
Fig. [5] while the stellar density profiles of the two satel- 
lites are well preserved at ~ 0.5 kpc, both profiles flatten 
on smaller scales due to interactions with the black hole. 
The typical radial scale on which we might expect the ef- 
fects of the black hole to appear is given by an equation 
analogous to M*(< tbh) = 2Mbh in Sec. 2.1.3, with the 
stellar mass referring to that of the satellite. This gives 
t'bh = 0.25 a for the satellites we are considering, i.e. 
~ 0.15 and 0.1 kpc for the la and 2a dense satellites, 
respectively, which indeed corresponds to the observed 
scale of the change in the density profile. We discuss the 
effects of black holes further in the next subsection. 

3.3. Dependence on presence of central black holes 

We have performed simulations without central black 
holes to quantify the effects of black holes in the host 
galaxies on the dynamics and disruption of infalling satel- 
lite galaxies. Fig. [6] compares results from runs with 
and without black holes for two orbits (left for radial 
orbit; right for v t /v r — 0.5; both with orbital energy 
^orb = l-5w c irc)- The black hole has a non-negligible ef- 
fect in the early evolution for the radial orbit because the 
satellite is brought to the host's black hole immediately. 
Both the central (dotted curve) and outer (dashed curve) 
regions are somewhat more disrupted in the case where 
the host has a black hole. The differences between the 
two runs become more pronounced as the merger pro- 
gresses: with successive pericentric passages the black 
hole exerts a stronger effect, widening the gap between 
the black and red curves on the left. By the time the 
merger is complete, the satellite stellar mass within its 
R e in the black hole run is only about half of that in the 
run without a black hole. The difference is even more 
pronounced at smaller radii; in fact, the black hole leads 
to a greater evacuation of mass within a than within R e , 
while the reverse is true for the run without a black hole. 

The runs with v t /v r — 0.5 (right panel) evolve differ- 
ently. The satellite does not pass near the black hole on 
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Fig. 5. — Evolution of the satellite's stellar density profile (top; logarithmic scale) and circular velocity profile (bottom; linear scale) for 
the same simulations shown in Fig. [4] Left: normal (SDSS) satellite. Center: la dense satellite. Right: 2a dense satellite. The profiles 
are computed at apocenter after each pericentric passage. Denser satellites are noticeably better preserved, both in their inner and outer 
structure. We do not plot the final circular velocity profile for the satellite because it is not meaningful when the satellite has merged with 
the host. Note the different radial scale on the upper and lower plots. 
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Fig. 6. — Comparison of the orbital decay and stellar mass loss of a satellite as it is accreted onto a host galaxy with (black) and without 
(red) a central black hole. Two initial orbits are shown: radial (left panel) and vt/v r = 0.5 (right panel). Black holes in host galaxies 
enhance the disruption efficiency for the satellite. 
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Fig. 7. — Evolution of the ratio of the satellite's stellar mass 
interior to the initial scale radius a sa t in runs with a black hole 
to the same quantity in runs without black holes. Two orbits are 
shown. Arrows mark first and second crossings of the black hole's 
sphere of influence tbh in the run with angular momentum. The 
ratio is seen to decrease noticeably between the crossings due to 
the tidal mass stripping by the black hole. The radial merger 
experiences crossings of reH at every pericentric passage, but the 
first 3 encounters are at high velocity and thus do not lead to 
significant tidal stripping. 

the first few pericentric passages, so the hole does not 
provide any shocking or stripping and there is no differ- 
ence between the runs with and without a black hole. 
Once the orbit has decayed enough to bring the satellite 
close to the center of the host bulge, however, the black 
hole starts to exert its effect. The final result is again 
that the satellite is more destroyed when the host has a 
black hole, although the effect is not as strong as for the 
radial collision. 

To quantify the amount of satellite destruction due to 
the host's black hole versus other tidal processes, we 
plot as a function of time in Fig. [7] the ratio of the 
satellite's stellar mass within its initial scale radius a, 
^W*,sat(< a , t), for the runs with a black hole relative 
to the same quantity in the runs without a black hole. 
For the radial merger (dashed line), this ratio decreases 
rapidly to about 25% after 0.4 Gyr due to the strong tidal 
effects exerted by the black hole on the accreting satellite. 
For the non-radial merger (solid line) , this ratio stays at 
essentially unity 4 until the satellite's orbit has decayed 
enough to cross the hole's sphere of influence (reH ~ 300 
pc) at ~ 0.5 Gyr. The first two crossings of tbh (marked 
by arrows) correspond quite well to the initial decrease 
in M* i8at (< a; BH) relative to M* iSat (< a;noBH). By 
the time the remnant has relaxed the run with a black 
hole has ~ 50% less mass within a sa t relative to the run 
without a black hole. This is less of a dramatic reduction 
than in the radial merger because the black hole inter- 
acts fewer times with the center of the satellite in the run 
with angular momentum. 

We have also considered the case where both the host 
and the satellite have central supermassive black holes. 
In this case, each black hole mass is set to 2 x 10~ 3 of 

4 There are oscillatory features from t ss 0.3 Gyr onward, but 
these simply reflect the slightly different orbital evolution of the 
runs with and without black holes; see the red vs. black solid 
curves in fig. \E\ 



its respective stellar bulge mass (see runs T5 & T6 in 
Table I). The dense satellite (with satellite-to-host mass 
ratio of 1:10) is used for this run; if there is additional dis- 
ruption of the satellite due to the presence of a black hole 
at its center, it will show up better in this case than in 
the lower-density one. We find that there is no substan- 
tive difference in the density profiles of either the host or 
satellite between the runs with one and two black holes. 
Our simulations, however, cannot follow the evolution of 
the system to the point at which three-body interactions 
among a supermassive black hole binary and stars start 
to scour out a core. The separation of the black holes 
at the hard binary stage is r^ard = G^/(A<7 2 ) ~ 1 pc 
for the parameters of our simulations; the highest force 
accuracy we use is e = 3 pc (run T6 in Table 1), and 
even this scale can be significantly affected b y two -body 
relaxation. Recent simulations by iMerrittl (|2006[ ) have 
studied the evolution of black hole binaries embedded in 
stellar spheroids with 7 profiles for a variety of black hole 
mass ratios and values of 7. The effect of a 1:10 binary 
(q = 0.1) on a 7 = 0.5 profile is not very pronounced 
(his Fig. 5), implying that the surface brightness profiles 
presented here would not be strongly affected by core 
scouring. 

3.4. Comparison to runs without dark matter 

Prior numerical studies on satellite accretion and 
their survival are focused on either purely dark mat- 
ter s u bhalos (e.g.. iHavashi et al.l 120031 : iKravtsov et al.l 
120041: iKazantzidis et all I2004J ) or purel y stellar sys- 
tems (e.g., the pioneering w ork o f iWhitd [198! 
iBalcells feQuinnl 119901: IWeinberd 119971: black holes are 
included in iHollev-Bockelmann fc Richstoriel 120001 and 
iMerritt fc Oruzll2001h . In this subsection we compare 
the results of simulations ignoring the effects of dark mat- 
ter to our full simulations (using stars, dark matter, and 
black holes) in order to quantify any potential differences. 
In general, simulations with only one primary component 
(stars or dark matter) are less expensive to perform than 
simulations with both components at a fixed force and 
mass resolution, so it is useful to assess whether including 
the effects of dark matter is necessary for understanding 
satellite evolution. 

As outlined in Sec. 12.11 there is in fact good cause 
for thinking that including both dark matter and stellar 
components will lead to non-trivial differences from simu- 
lations with only dark matter halos or stellar bulges. The 
results of the previous sections indicate that one impor- 
tant parameter in determining whether or not a satellite 
is strongly disrupted is the number of pericentric pas- 
sages it experiences before merging with the host bulge. 
Varying the depth of the host potential well is likely to 
change the number of orbits a satellite undergoes. Since 
orbital deceleration is governed by dynamical friction 
(equation [2|) , which is proportional to M sat (i) / 'Mhost(^) , 
and tidal stripping is likely to reduce M sat by removing 
dark matter from outside the tidal radius, we expect the 
simulations with dark matter to have a slower orbital 
evolution and to go through more pericentric passages 
before merging than the equivalent runs without dark 
matter. 

To quantify the effect of dark matter, we compare in 
Fig. [5] the results from a 3-component run (with dark 
matter, stars, and black hole) and a 2-component run 
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Fig. 8. — Comparison of radial satellite accretion with (left) and without (right) dark matter halos. The run including dark matter (in 
both satellite and host) has ti or t, = 0.75t> c i rc while the run without dark matter has v OI ^ = v c i rc . Based purely on the results presented to 
this point, we expect the run with the higher ratio of f or b/t>circ to be more disrupted. The reverse is true, however, because the dark matter 
significantly deepens the host galaxy's gravitational potential, leading to more pericentric passages for the satellite with dark matter and 
hence more severe mass losses. 

the r~ 1 cusp. For example, our standard density satel- 
lite (M* = 3.3 x 10 10 M Q , a = 0.76 kpc, 7*, s = 1) looks 
quite similar to a satellite with the same stellar mass but 
7* B = 1.5 and a = 1.35 kpc: the densities differ by a 
maximum of 30% from 30 pc (our standard force soft- 
ening) to 5 kpc. We include analyses of the runs with 
steeper satellite stellar density in the sections that follow 
below. 



without the dark matter halo (run LI and Tl in Table 1). 
(The host and satellite either both have dark matter ha- 
los or both without.) The satellite for both runs is on a 
radial orbit with initial orbital velocity u or b = 0.75i; c i rc 
(at 4R e ). Note that the circular velocity v c i YC is not equal 
in the two cases, as the dark matter contributes to the 
circular velocity. Fig.[8]shows that the satellite in the run 
with dark matter (left panel) experiences more pericen- 
tric passages before the final merger than the run with- 
out dark matter (right panel). Even though the satel- 
lites' stellar masses interior to R e and a are quite similar 
between the two runs after the first pericentric passage, 
the additional ~ 3 passages in the run with dark matter 
cause extra disruption. As a result, only 20-30% of the 
satellite's stellar mass remains within its initial R e and a 
after 0.2 Gyr in the 3-component run (dashed and dotted 
line in left panel) compared with the 70-80% in the run 
that ignored dark matter. 

3.5. Other factors: mass ratio, inner density profile 

We have tested a number of other parameters that can 
potentially affect the outcome of our satellite accretion 
simulations. Primary among these are the satellite-to- 
host mass ratio and the inner logarithmic slope 7* !S of 
the satellite's density profile. 

All the results presented thus far have assumed a 
satellite-to-host mass ratio of 1:10. We have performed a 
two-component test simulation (run T2 in Table 1) with 
a lower mass ratio of 1:33 and found qualitatively similar 
results. These lighter satellites lose their orbital energy 
on a longer time scale, however: it takes a factor of ~1.5- 
2.5 longer for their stellar bulges to merge with and have 
observational impact on the host's bulge. In a cosmolog- 
ical setting, we therefore expect the satellites with < 5% 
of the host mass to have less direct dynamical effect on 
the host than the more massive satellites despite their 
higher number density. 

We have also performed several simulations (runs la- 
beled "s" in Table 1) using a steeper p oc r~ 15 cusp for 
the initial stellar density of the satellite. As expected, 
these satellites suffer less tidal disruption than the r _1 
satellites studied thus far, but the net result is somewhat 
degenerate with decreasing R e of the satellite while fixing 



4. OBSERVABLE CONSEQUENCES OF SATELLITE 
ACCRETION 

4.1. Surface brightness profiles 

As we have shown in Sec El the fate of a sinking satel- 
lite depends sensitively on a number of parameters in- 
cluding its orbital energy and angular momentum, its 
mass concentration, and the presence of black holes and 
dark matter halos. Even one dense satellite galaxy that 
finds its way to the center of a massive elliptical galaxy 
without being significantly disrupted can lead to a cen- 
tral density enhancement that is occasionally observed, 
so surface brightness profiles can serve as a useful diag- 
nostic for structure formation models. 

In Fig. [9] we illustrate how a single sinking satellite 
on a variety of orbits and with a range of stellar densi- 
ties can modify the shallow surface brightness profile of 
a massive elliptical galaxy described in Sec 2. Fig. fTOl is 
a busier summary plot for all of our 3-componcnt simu- 
lations listed in Table 1 (including those in Fig. [5]). Each 
curve is computed using the mean of twenty random pro- 
jections of the merger remnant. In Fig. [9l the upper left 
panels compare a relatively bound satellite orbit at fixed 
energy (Worb/fcirc = 0.75) with varying angular momen- 
tum; the upper right panels show a more energetic or- 
bit (worb/w circ = 1-5) with varying angular momentum; 
the lower left panels compare runs with identical orbital 
parameters (w rb/w C irc = 1-5, Vt/v r — 0.5) but varying 
initial satellite densities; and the lower right panels com- 
pare energetic (vorb/wcirc = 1-5) mergers with steeper 
satellite stellar density (7*, s = 1.5) and varying angular 
momentum. The top panels show the absolute surface 
brightness profile S(i?) of the stellar component in the 
host galaxy initially (thin solid lines) and after satellite 
accretion (other lines; computed from all stellar particles 



Satellite Accretion Onto Massive Galaxies 



11 




0.1 1.0 0.1 1.0 

R (kpc) R (kpc) 




0.1 1.0 0.1 1.0 

R (kpc) R (kpc) 

Fig. 9. — Upper of each set of panels: Surface brightness profiles S(-R) of the host galaxy initially (thin solid curves) and after accretion 
of a satellite (computed from all stellar particles) for selected three-component runs. Lower of each set of panels: Same as top panels 
except it shows the corresponding fractional change between the initial and final 2(iJ) to accentuate the evolution of the surface brightness. 
Upper left: runs with varying orbital angular momentum at fixed orbital energy of ^ or b/«circ = 0.75 and 7*, s = 1.0. The final central 
surface brightness is larger than the initial for all three runs but the amount of increase decreases with decreasing angular momentum. 
Upper right: runs with varying orbital angular momentum at fixed orbital energy of forb/^circ = 1-50 and 7* jS = 1.0. The increased orbital 
velocity leads to less deposition of satellite mass at the center of the remnant than in the v OI ^/v c i IC = 0.75 case, and the radial collision 
results in a noticeably reduced inner central surface brightness. Lower left: runs with fixed orbital angular momentum (vt/v r = 0.5) and 
fixed orbital energy (v orh /v c - lIC = 1.5) with 7* jS = 1.0 and varying satellite stellar density (i.e. varying R e )- The standard density satellite 
results in a reduced central surface brightness while the more compact satellites have dense centers that survive to the center of the host 
bulge, adding a large amount of mass. Lower right: runs with varying orbital angular momentum at fixed orbital energy of v OI ^/v c - lIC = 1.5 
and 7*, s = 1.5. The final central surface brightness decreases for the radial run but increases for both runs with angular momentum. 
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Fig. 10. — Fractional change in the central surface brightness 
profile of a host galaxy after accretion of a satellite from all 3- 
component simulations listed in Table 1. The profiles show a wide 
spectrum of behavior on small radial scales, ranging from reduction 
in the central surface brightness for energetic and eccentric orbits, 
to sharp increases when the accreted satellite is relatively compact 
and dense initially. All the final surface brightness profiles have 
similar behavior at larger radii (> 2 kpc), regardless of the central 
profile. 

in both the host and satellite); the bottom panels show 
the fractional changes between the initial and final 

The general trend in the upper two sets of panels in 
Fig. [9] and in Fig. [10] is that satellites on more energetic 
or more eccentric (i.e. lower angular momentum) orbits 
are more destroyed due to the stronger tidal effects. This 
trend is consistent with the results in Sec. 13. H and Figs. [2] 
and [3] While Figs. [2] and [3] illustrate the orbital decay 
and mass loss of the satellite, Figs. [9] and \W\ show the 
net effect of the sinking satellite on the inner stellar dis- 
tribution of the host galaxy. 

A particularly noteworthy feature is seen in several 
cases in Fig. [S] and in Fig.[TUJ the central surface bright- 
ness of a core elliptical galaxy is actually reduced when a 
satellite is accreted on a relatively energetic (t> or b ^ w c j rc ) 
and eccentric (vt < v r ) orbit. Satellite accretion in this 
particular orbital parameter space is likely to occur infre- 
quently, but it provides a potential mechanism for pro- 
ducing a minimum in the central surface brig htness. Min- 
ima are see n in ~ 5% of the core galaxies (jLauer et al.l 
2002, 2005), although their physical scale tends to be at 
< 100 pc. A broader search in parameter space would be 
needed to determine whether the observed features can 
be more closely reproduced by specific combinations of 
parameters in this satellite accretion scenario. 

The somewhat counter-intuitive reduction in central 
E(i?) via satellite accretion is the result of two processes. 
First, the accreting satellite clearly has to be largely de- 
stroyed so as to bring in little stellar mass to the host's 
center. As Figs. [2] and [3] showed, energetic and eccentric 
orbits are effective means of destroying satellites. Sec- 
ond, as a satellite sinks via dynamical friction toward 
the host's center, energy is transferred from the satel- 
lite to the host bulge, leading to gravitational heating 
that helps reduce the central stellar density in the host. 



Such satellite-host dynamical interaction and its effect on 
the bulge density and surface brightness profiles is there- 
fore similar to that for dark matter halo density profiles, 
where dark matter subhalo accretion tends to heat the 
central cusp, reducing the central density, while depo- 
sition of the stripped mass and merging of the subhalo 
tends add to the central density (|Ma &: Bovlan-Kolchinl 
l200l . In both cases (satellite accretion and dark mat- 
ter subhalos) it is the competition between gravitational 
heating and mass deposition that determines whether the 
final profile is more dense or more diffuse than the initial 
profile. 

The lower left panels of Fig. [9] show the strong effect 
of the initial satellite density on the host's central S(i?). 
The satellite's orbital decay, mass loss, and the evolution 
of its interior structure from the same three simulations 
were shown in Figs. 2] and [5l While the satellite on the 
SDSS R e — M± relation (thick solid lines) is disrupted and 
leads to a ~ 10% overall reduction in S(i?) in the inner 
~ 500 pc (see paragraph above), the denser satellites 
(dashed and dotted lines) lead to a dramatic increase in 
the central surface brightness. Similar results are seen for 
satellites with 7* !S = 1.5 rather than 1.0 with substan- 
tial angular momentum (lower right panels). Accretion 
of dense satellite - either with a steeper central stellar 
density profile or with a smaller R e - is therefore likely 
to add a dense remnant to the center of a core galaxy 
even when the tidal effects from the black holes are in- 
cluded. A small number of surface brightn ess profiles 
from the galaxy sample in lLauer et al.l (|2005f ) somewhat 
resemble the profiles of the 1 — 2a dense satellites (e.g., 
NGC 3384). We suggest that satellite accretion offer a 
pla usible explanation for the shapes of those profiles (see 
also lKormendviri984f ) . 

4.2. Colors 

Both the average color and the color gradients in el- 
liptical galaxies contain vital information about their 
star formation and assembly histories. For instance, the 
slope and scatter along the red sequence in the color- 
magnitude relation (CMR) is a useful diagnostic of how 
pre-existing early-type galaxies build up larger ellipti- 
cal galaxies, as this merging would tend to average the 
color s and thus wash out any correlations ([Bower et alJ 
1998). Monolithic collapse models generally predict steep 
color gradient s originating from strong ra dial metallic- 
ity gradients (|LarsorJll9"74T; ICarlberdll984h . The case is 
not as clear-cut in hierarchical models, as the effects of 
metallicity gradients established by gas-rich mergers can 
be diluted by violent relaxation and mixing in gas-poor 
mergers. In general, the distribution of a satellite's stars 
should affect the local color of a host in a way that de- 
pends on the satellite orbit and density. If the satellite 
sinks to the host's center without substantial mass loss, 
for example, the resulting color on small scales will be 
close to the satellite's initial color. If, on the other hand, 
the satellite is disrupted outside the host's core, the rem- 
nant's central color will be essentially identical to that of 
the host initially. 

We investigate the changes in galaxy colors due to 
satellite accretion by assigning initial colors to the host 
and satellite and then computing the remnant's color un- 
der the assumption of no change in color due to stel- 
lar population evolution. This method isolates the ef- 
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Fig. 11. — 01 {g — r) color profiles for mergers with dark matter and a black hole, each averaged over 20 random projections. The host is 
assumed to have a color of 1.0 initially, while the satellite is given a color of 0.9 in each case. We consider two cases: no color gradient in 
either galaxy (red curves) or a color gradient of -0.1 mag/dcx in both satellite and host (black curves). The panels are arranged as in Fig. [9] 
Upper left: v OI \, = 0.75v c i IC with varying vt/v r . Upper right: it OI b = 1.5 v c irc with varying vt/v r . Lower left: v OI ^, = 1.5u c i ro , Vt/vr = 0.5 
with varying satellite density. Accretion of a compact satellite leads to a relatively bluer central remnant, while standard density satellites 
cause less change regardless of orbit. Lower right: i> or b = 1.5u c irc an d 7*,s = 1.5 with varying vt/v r . 
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Fig. 12. — Two-dimensional color map of the host bulge, viewed in the orbital plane, after accretion of a satellite with different 
initial densities: standard (left), la denser (middle), and 2a denser (right). The satellite orbit is identical for all three cases: vt/v r = 
0.5, forb/^circ = neither galaxy is given an initial color gradient. Each panel is (6.66 kpc) 2 ; each pixel corresponds to 0.1 kpc. The 
cores of the 1 and 2a dense satellites survive to leave a visible core bluer by as 0.05 — 0.08 mag. Close inspection of the left panel shows 
that although the satellite with standard density does not survive to the center, it leaves faint plumes of stripped stars (light red) at very 
low color contrast (as 0.01 mag bluer). 
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fects of mixing in the mergers. We assign the initial 
galaxies 01 (5 — r) colors 5 of 0.9 for the satellite and 1.0 
for the host. The observed CMR slope from SDSS in 
01 {g - r) is -0.022 mag/dex ()Hogg et alj |2004), meaning 
we have slightly accentuated the expected color differ- 
ence between the host and satellite given their masses. 
When we consider galaxies with initial color gradients, 
this color corresponds to the color measured within the 
effective radius. 

Fig. [TT] shows the circularly-averaged color profiles for 
the same sets of simulations shown in Fig. [9] as in Fig. [9j 
each curve is computed using the mean of twenty ran- 
dom projections of the remnant. The red curves assume 
no initial color gradient for either galaxy and show that 
the properties of the satellite and its orbit make a sig- 
nificant difference in the color of the inner ~ 1 kpc of 
the remnant. The dense satellites (lower left panel) and 
7*. s = 1.5 satellites (lower right panel) lead to a bluer 
center, a result of their undisrupted cores settling to the 
center of the host. The runs with standard density and 
varying orbital properties show less pronounced central 
differences. For all simulations, the 01 (g — r) color at 
large radii is about 0.99, so A color w 0.01 mag. 

The black curves in Fig. [TJJ assume each galaxy ini- 
tially had a color gradient of —0.1 mag/dex and show 
how these gradients evolve as a result of the merger. In 
the runs with lower w or b (upper left panel), the gradi- 
ent is diluted somewhat in the center but only the run 
with the highest angular momentum causes a central in- 
version. None of the runs at higher v OI b and 7* jS = 1.0 
(upper right panel) have a central color inversion, but 
the color profile does flatten toward the center. The 
runs with higher stellar density or steeper inner stellar 
density slope show a marked change in central color at 
approximately R e /5 (as in the case without initial gra- 
dients). The general trend we see is to weaken the cen- 
tral color gradients while affecting the outer regions less 
strongly. T his wea k ening is in agreement with the trend 
reported in I White! (|1980l ) . Luminous elliptical galaxies 
do tend to have slig htly weaker gradients in meta llicity 
(| Car olio et al.lfl993l ) and color (jLauer et al.l l2005h than 
faint ellipticals galaxies, perhaps as a result of this mix- 
ing process. 

To quantify the color signatures of satellite accretion 
further, we have also examined two-dimensional color 
maps of the hosts to look for features that are not evi- 
dent in the 1-dimensional color profiles shown in Fig. 1111 
A set of such maps (viewed in the orbital plane) is shown 
in Fig. [T2l comparing the three runs with different initial 
satellite densities. The initial colors are the same as for 
Fig. [TTJ The 2a dense satellite (right panel) sinks to the 
center without significant disruption, leaving an observ- 
able core of w 0.5 kpc in the remnant. Even the la dense 
satellite (center panel) leaves an observable color signa- 
ture in the core, although it is both smaller and and less 
concentrated than the 2a dense core. By contrast, the 
standard density satellite is largely destroyed and leaves 
no core in the remnant (left panel of Fig. fT2|) . as was 
previously seen in other ways (e.g., Figs 0] and O . Close 
inspection, however, shows faint plumes of bluer features 

5 The exponent 0.1 in the (g-r) color indicates that the SDSS 
galaxy magnitudes are K- corrected to z=0.1 rather than z=0.0 - 
sec, e.g. IHoeg et ail l[2u0l ). 



at low color differences (« 0.01 — 0.02 mag) from stripped 
stars of the torn-up satellite. The axis of the plume cor- 
relates well with the orbit of the satellite over its final 
several pericentric passages (which are quite radial). 

4.3. Additional observable properties: intra-cluster light 
and kinematic structure 

Many massive galaxies at or near the centers of clus- 
ters are cD galaxies with extended stellar envelopes, 
often extending out to hundreds of kpc. Satellite ac- 
cretion onto cluster halos (and subsequent satellite dis- 
ruption) contributes, p erhaps significantly, to this intra- 
cluster light (ICL; e g lMerrittlfl984l: iMoore et all 119961: 
[ Gregg fc West! 119981 : iGnedinl 120031 : IZibetti et all 120051 : 
iMihos et al.ll2005T ). Our initial setup for the galaxies is 
biased against producing ICL in the sense that we trun- 
cate the stellar bulge, removing the stars that are at large 
radii (and thus are least bound). This truncation elimi- 
nates the stars that would be most easily removed from 
the satellite and deposited at large radii relative to the 
center of the remnant. Nevertheless, it is not inconceiv- 
able that a fraction of stars in our mergers is deposited 
out to large radii, contributing to an ICL component. 

The remnants in all of our simulations show excess stel- 
lar mass at large radii (> 3R e for the initial host profile), 
as is visible in the lower panels of Fig. [9] and Fig. [10] Al- 
though the excess is relatively small - of order 15-25% - 
it shows that the star particles do get deposited signif- 
icantly outside of the host's effective radius, indicating 
that tidal stripping of merging satellites can contribute 
to ICL. For cases with w or b = 1.5t> C i rc , the stellar mass 
that belonged to the satellite tends to be distributed with 
a larger characteristic radius than the mass that belonged 
to the host. Simulations of repeated satellite accretion in 
a cluster environment would likely yield a significant ICL 
component originating from the tidal stripping of satel- 
lites. Analyses of cosmological si mulations have shown 
this process to yield some ICL (|S ommer-Larsen et al.l 
2005? iWillman et all l200l iRudick et all 120061) : future 
higher-resolution simulations should shed more light on 
the origin of the ICL. 

Many elliptical galaxies are observed to have kinemat- 
ically distinct substructure at their centers (kinemati- 
cally decoupled cores, KDCs); often these cores counter - 
rotate relative to the galaxy (|de Zeeuw fc Franxlll991h . 
A recent survey by the SAURON team of 48 nearby 
(cz < 3000 kms" 1 ) E/S0 galaxies revealed a wide range 
of kinematic substructure, from ~ 100 pc young cores 
in fa st rotators to kpc-scal e, old KDCs in slow rota- 
tors (|McDermid et al. 2006). We have tested for such 
substructure by making velocity maps of our remnants. 
None of them exhibit a KDC; the velocity fields in all 
cases are consistent with zero net rotation and no kine- 
matic subcomponents. Intriguingly, the three most lu- 
minous galaxies in the SAURON sample, including M87, 
do not have KDCs. Since this type of galaxy is similar 
to our initial host, our simulations are in this sense in 
good agreement with the SAURON observations. 

A broader question to be explored in future work is 
the origin and fragility of KDCs in slow-rotating galax- 
ies. Gas physics co u ld pla y a role in cre a ting KDCs: 
Hernauist fc Barnesl (jl991l ). Uesseit et ail (|2006l ). and 
Cox et al. 1 (12006ft tiave demonstrated that simulations of 
binary disk galaxy mergers with gas can lead to kine- 
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matic subst ructure in the remnants . In dissipationless 
simulations, iBalcells fc Quinrjj |1990|) have shown that a 
dense satellite sinking to the center of a more diffuse 
host can create a counter-rotating core if a large initial 
rotation of V/<jq — 0.75 for the satellite and 0.4 for the 
host are used and the merging orbits are retrograde. It 
is somewhat difficult to discern from the SAURON re- 
sults what fraction of slow-rotators should have KDCs 
given the sample size (12 slow- rotators) . If their results 
are representative, however, then it is possible that any 
KDC would be erased if massive ellipticals are built by 
multiple generations of merging. 

One satellite accretion event with mass ratio 1:10 does 
not strongly affect the radius or velocity dispersions of 
the remnant nor does it increase M* substantially, mean- 
ing scaling relations such as the fundamental plane or its 
projections do not change notably. Remnants of major 
mergers (mass ratio 1:1-1:3), on the other hand, can have 
very different properties dependin g on energy or angular 
mom entum of the merger orbit ( Bovlan-Kolchin et alJ 
2006). Checking the evolution of sizes, stellar masses, 
and velocity dispersion over a sequence of merger events 
would therefore be a good test of the evolution of scaling 
laws such as the Faber- Jackson and Mbh — c e relations. 

5. CONCLUSIONS AND DISCUSSION 

In this paper we have used dissipationless simulations 
to examine in detail and over a wide range of parame- 
ter space the case-by-case outcome of a single satellite 
elliptical galaxy accreting onto a core elliptical galaxy. 
We have taken special care to explore cosmologically rel- 
evant orbital parameters and to set up realistic galaxy 
models that include all three relevant dynamical com- 
ponents: dark matter halos, stellar bulges, and central 
massive black holes. The main results in this paper are 
as follows: 

(i) Satellites on more energetic or more eccentric (lower 
angular momentum) orbits are destroyed more efficiently 
due to the stronger tidal effects from the host spheroid 
and the central black hole (see Figs. 2, 3, 9-12). 

(ii) Denser satellite galaxies - those with either a 
smaller R e or a steeper inner slope - survive significantly 
better (see Figs. HI 5, 9-12). For instance, the standard 
density satellite in Fig. [4] loses more mass by its 3rd peri- 
centric passage than the total mass loss in the 2a dense 
satellite during its entire orbital evolution. 

(iii) A supermassive black hole at the center of the host 
galaxy significantly increases the efficiency with which 
an accreting satellite is destroyed (see Figs. 6, 7, 10). 
Without a central black hole, luminous elliptical galax- 
ies would have more difficulty in protecting their shal- 
low inner profile (however formed initially) from steep- 
ening due to mass deposit through satellite accretion, in 
agreement with previous 2-component numerical work 
(iHollev-Bockelmann fc Richstonel l2000t iMerritt k. Cnizl 
1200 ID . 



(iv) Ignoring dark matter halos in the galaxy models 



causes significant differences primarily because dark mat- 
ter halos lead to higher orbital velocities for a fixed ratio 
of Worb/wcirc (see Fig. 8). Satellites with increased or- 
bital velocity complete more orbits before merging with 
the host, resulting in greater destruction of the satellite 
galaxy. 

(v) Our overall conclusion is that the accretion of a 
single satellite elliptical galaxy (1/30 to 1/10 mass of the 
host) can result in a broad range of changes, in both 
signs, in the surface brightness profile X(-R) and color of 
the central part of an elliptical galaxy, as summarized in 
Figs. [101-12. 

More specifically, Figs. [T0lT2 show that when a black 
hole is included in the host galaxy and when the satellite 
and host are both on the SDSS stellar size-mass relation, 
the fractional change in S(i?) in the inner ~ 500 pc of the 
host galaxy upon satellite accretion ranges from ~ —20% 
to ~ +40% When more compact satellites are consid- 
ered in accordance to the large scatter about the SDSS 
size-mass relation, the inner S(i?) of the host galaxy can 
increase by more than 100%. Since accretion of a sin- 
gle satellite adds little mass and luminosity to the host, 
we suggest that a "core" galaxy at M v ~ —22 can po- 
tentially transition into a galaxy with power-law char- 
acteristics at a comparable luminosity by accretion of 
a relatively compact and massive (~ 1:10) satellite. A 
core-in-core structure would prov ide particularly s trong 
evidence for this process (see also iKormen dvl fT98l . 

The strong dependence of satellite mass loss on the 
satellite's density relative to the host and the inclu- 
sion of black holes has broader implications, for exam- 
ple, in interpreting the connection between dark mat- 
ter subhalos and galaxies in cosmological simulations 
that model only dark matter. Dark matter subhalos 
in these simulations are often stripped to the point 
of not being distinguishable from the host dark mat- 
ter h alo. In semi-analytic models of g alaxy formation 
(e.g. iGao et al.|[200l I Wang et al.ll2006l i. these subhalos 
are tracked as "galaxies" by either following the most 
bound particle of the subhalo or extrapolating their po- 
sitions and velocities using information from previous 
timcstcps. The results from our three-component (dark 
matter +stars+black hole) non-cosmological simulations 
suggest that the properties of the observable stellar com- 
ponent in the simulations are very sensitive to their as- 
sumed initial properties, and in particular, to the initial 
stellar density assigned to the galaxies and the inclusion 
of black holes in the host galaxies. Future work that 
incorporates the results presented here into a full cosmo- 
logical context should be valuable. 
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APPENDIX 

PROPERTIES OF TWO-COMPONENT 7 MODELS 
Untruncated Models 

iDehnenl (|1993l Appendix B) gives several useful formulae for numerically computing quantities such as surface 
brightness profiles and projected velocity dispersions for standard 7 models. His integral transforms can be readily 
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Fig. 13. — Left: Density profiles of 7-profiles used in this paper. Solid curves show the untruncated profiles (7* h = 0.5, M = 
3.33 X 10 11 Mq, a = 2.33 kpc for the bulge, 7dm = 1, M = 3.33 X 10 12 M , a = 33.3kpc for the halo). The dotted lines show truncated 
profiles for the halo (from left to right, a,' DM = 66.7, 133, 333 kpc; 133 kpc is used for the host halos in the simulations), while the 
dot-dot-dot-dash lines show truncated profiles for the bulge (from left to right, a'^ = 33.3, 83.3, 167 kpc; 167 kpc is used for the host bulges 
in the simulations). Note that as a' — > a (leftmost dotted curve), the inner density amplitude drops but the slope remains the same. Right: 
Stability test of a standard satellite for our simulations. Plotted are the initial density profiles (solid curves) for the stellar bulge and dark 
matter halo, as well as the evolved profiles at t = 1 Gyr (dashed curves). The force softening for this run is e = 33 pc. The model is stable 
on all scales larger than ~ 2.5 e. 



extended to both full and truncated two-component 7 models, as we show here. We first consider the untruncated 
models composed of three components (dark matter, stars, and a central black hole), with /?dm = /?(7dm, A/dm, «dm; r) 
and p± = p(7*, M*, a*; r); we set Newton's constant G to unity in the calculations that follow. 

Since the spatial structure of each component is unaffected by the other component, the formulae in iDehnenl (|1993l ) 
for surface density and cumulative surface density are valid for each component of the untruncated two-component 
models. Computing the projected velocity dispersion is somewhat more complicated, as the velocity structure of a 
system is determined by the total gravitational potential. 

The projected stellar velocity dispersion <j prk (R) as a function of projected radius R can be computed from the Jeans 
equation: 

2 f°° r 

"'■ JR) = K(R)L "•"' <r) (r'-R'yi' dr < A1 > 

(c.f. iTremaine et~aTlll994L eqns. 27-29), where c 2 (r) is the stellar radial velocity dispersion. The total mass profile can 
be split into three parts, M tot (r) — M*(r) + A/dmM + A/bhMj meaning the integral can be split into three pieces: 
°p,*( r ) = s p,*( r ) + Sp DM (r) + Sp BH (r). The first term is the projected dispersion for the one-component case (e.g., 
Dehnen's eqn. 18), though a singularity in the integrand means it can be much mor e easily computed usin g his eqn. B2. 
The third term can be mapped into an integral similar to that of the first term (|Tremaine et al.lll994 eqns. 49-51) 
and then computed using Dehnen's eqn. B2. This leaves only one piece to be computed: 

4dmM = ^ [ ElM^H^ _ R ^ dr . (A3 ) 

This term is difficult to compute numerically because of the limits of integration. Following Dehnen, we make the 
substitution 

t 2 = ^— {8 + 1) -8, (A4) 
r + a* 



with s — R/a+, and take a* = 1 = M*. Eqn. IA3I then becomes 

s 2 (J?)- MdmM* (5 : , , 



x / (1 - t 2 f{t 2 + fl )i-T*-TD M [ t 2 (1 _ aDM ) + s + aDM ]yoM-3 t 2 ^ 2s + f2 ^ _ s) dt ^ (A5) 



and numerical integration is straightforward. 
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Truncated Models 

As described in the main text, we truncate the density profiles of our models by subtracting off a second 7 profile: 
assuming the desired density profile is Pi(j, M, a; r), we take the truncated density profile to be 



Pi(r) = Pi(7, M, a- r) - p t {i , M' , a'; r) . 



(A6) 



In order to preserve the same inner density structure, we require that 7' = 7, M' < M, and a' 3> a. By varying M' 
and a', we can control the truncation. The maximum radius for one component in these models is given by 



where 



A = 



M 
M 7 



A-l 

''1-4A 



3-vH V(4-7) 



(A7) 



(A8) 



The left panel of Fig. [T3] shows the profiles, both untruncated and truncated, for the dark matter halos and bulges 
used in several of the simulations presented in this paper. The numerical stability of one of the satellite models is 
demonstrated in the right panel of Fig. [T2J there is no evolution in the density profile for scales > 2.5 e. The host 
models are equally stable, both with and without black holes. 

Since the truncated models do not have infinite extent, many of the integration limits change. For example, the 
velocity dispersion has f — * / J , max , where r max is given by eqn. IA7I Projected quantities therefore have an upper 
limit of r max rather than 00, or equivalently, £ max = t(^max) rather than t max = 1. The calculations are also changed 
by the addition of the second profile, which essentially serves as two additional mass components. For example, now 
Mtot — ► AI* — M*i + A/dm — ^Idw + Mbh and p*(r) — > [p*(r) — P*'(r)] (and similarly for pbm)- Velocity dispersions 
can be computed substituting these relations into the appropriate equations given above. 

The surface brightness for a truncated model is similarly given by 



p(r) 



rdr 



( r 2_ #2)1/2 

so introducing a truncation for the bulge in the form of p*>(r) = p(7*' 5 M*i, av; r) is simply a linear operation and 



(A9) 



= E*(i?) -Ev(iZ). 



(A10) 
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